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Emitted radiation can be reprocessed in curved space- 
times, due to the breakdown of the Huyghens principle. A 
maximization procedure for the energy diffusion allows one to 
obtain wave packets (gravitational and electromagnetic) that 
are particularly strongly backscattered. Examples are shown 
with the backscattered part exceeding by one order remnants 
of initial signals. A robust ringing can be observed, with am- 
plitudes exceeding leftovers of the main radiation pulse. The 
analysis of the obtained results allows one to set demands on 
some parameters in the numerical description of a realistic 
process of the collapse of two black holes. 



I. INTRODUCTION 

It is known essentially since the times of Hadamard 
[1] that curved spacetimes can affect the propagation of 
waves. The breakdown of the Huyghens principle [1] (or 
the backscatter, a name adopted by the general relativity 
community after de Witt and Brehme [2]) can influence 
both the energy and the energy flux of a wave signal. 
The backscatter can leave its imprints on the frequency 
spectrum and can affect the transmission time. The man- 
ifestations of this effect are the so-called tails and, most 
impressively, the quasinormal modes (QNM thereafter). 
The literature on the backscattering and related phenom- 
ena is quite extensive - see ( [3] ~ [23]) and numerous 
references therein. 

The QNM's have some features of the scattering- type 
solutions and they have been studied in the context of 
general relativity for more than three decades [6] . Many 
of their characteristics are well known for black holes 
[14], for instance their (complex) frequency spectrum. 
An observer located at a fixed space position would find 
that QNM's oscillate with amplitudes decreasing expo- 
nentially in time. The oscillation periods and the damp- 
ing exponents are the real and the imaginary parts of 
a frequency, respectively. They depend only on a few 
global characteristics of black holes - their asymptotic 
mass, angular momentum and/or global electric charge. 
Therefore their identification in an observed wave spec- 
trum would unambiguously identify a black hole (and in 
fact provide an argument, closest to the direct observa- 
tion, in favour of the existence of black holes) . Extensive 
reviews are presented in [20] and [21]. Tail terms have 
been studied in 1970's beginning from Price [8] but in- 
terest in them has again revived recently [23] . 



The spectra of QNM's and the decay exponents of the 
tails are universal, independent of initial data, but the 
very existence of QNM's and their amplitudes (as well of 
the tails) do depend on initial wave conditions. The main 
aim of this paper is to show the strongest imprints of 
the spacetime curvature that are present in the 
form of QNM's in a propagating wave. The im- 
plementation of this task requires the separation of the 
genuine geometric effects from those being built into ini- 
tial data - notice that even in the Minkowski spacetime 
one can easily form a QNM-like structure by produc- 
ing suitable initial data. The simplest possibility is to 
consider the purely backscattered part of the initial ra- 
diation, which is absent in the Minkowski geometry but 
which always exists in a curved spacetime. 

It would be meaningless to try to accomplish our aim 
by the method of " trial and error" ~ by selecting at ran- 
dom various initial wave configurations from the ocean 
of all possible data. Rather one should focus on "ex- 
tremal" in some sense initial data that can generate, in 
the first instance, "extremal" asymptotic templates, but 
also can set bounds on some parameters used in the nu- 
merical descriptions. In the present paper we follow the 
second strategy, using as a guiding principle the idea of 
extremizing the so-called diffusion parameter [24] and ad- 
dressing following issues. First, we estimate the maximal 
strength of the backscatter. The corresponding profiles 
of initial wave packets are found to favour vigorous ring- 
ing and/or strong deformation of initial signals. Second, 
and in relation with the former point, we obtain infor- 
mation on the process of taking waveforms from specific 
properties of the backscattered radiation. The order of 
the rest of this paper is as follows. Sec. 2 provides ba- 
sic information on the wave equations. Sec. 3 describes 
in detail the procedure of maximizing the diffusion pa- 
rameter and shows exemplary initial data for the wave 
evolution. Sec. 4 reviews some representative examples 
of wave templates. In Sec. 5 we again review those fea- 
tures of the numerical examples that could be useful for 
the numerical relativists dealing with the full nonlinear 
description of the collapse of two black holes. Sec. 6 
summarizes main conclusions. 
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II. BASIC DEFINITIONS AND CONCEPTS 



A. Equations 



The spacetime geometry is defined by the hne element 



ds-" = -rjRdr + + R-'dfl'', 

Vr 



(1) 



where t is a time coordinate, R is the radial areal coordi- 
nate, 7?i^ = 1 - ^ and dfl"^ = dO'^ + sin^ Odc/)'^ is the line 
element on the unit sphere, < (j) < 2tt and < 9 < tt. 
Throughout this paper the Newtonian constant G and 
the velocity of light c are put equal to 1. 

We will study the propagation of polar and axial 
modes of the quadrupole gravitational waves (GW there- 
after) and the dipole electromagnetic waves (EW) in the 
Schwarzschild background. The evolution equation has 
the form 



(2) 



Here r* = R + 2m\n^-^ — ij is the tortoise coordinate 
while the potential term reads: for the polar GW 



i?2 ^ '"2i?4(l + |m)2' 



for the axial GW 



and for the dipole EW 



V{R) = 2^. 



(3) 



(4) 



(5) 



The evolution equations corresponding to the first two 
potentials are called the Zerilli equation [25] and the 
Regge- Wheeler equation [26], respectively. 



B. Conserved energy 



The equation (2) possesses a conserved energy, 



E{R,t)= J ^(^{dt^f + {dr^^f + V^^y, 



(6) 



that is, the rate of change of in a fixed volume equals 
to the total flow through the boundary ( [27] and [28]). 
This agrees (up to a constant factor) with the energy de- 
duced from the stress-energy tensor for the EW. Eq. (6) 
represents a mathematically useful quantity in the case of 
the gravitational waves, with the density being asymp- 
totically proportional to the density in the quadrupole 



formula. In cither case, the energy conservation becomes 
important in our forthcoming construction. 

Assume that initial data ^ and dt'^ vanish inside a 
sphere having a radius a > 2m. From the conservation 
law one easily finds that the amount of the energy that 
reaches a distant observer is equal to 



where 



Ea{oo) = E{a,t = 0) - dEa, 



6Ea = dt[{dt^ + a^.*)2 + 



(7) 



(8) 



( [27] and [28]). The integration in (8) is done along the 

outgoing null cone that starts from a at t = 0. In the 
Minkowski spacetime (put formally to = in Eq. (2)) all 
of an initially outgoing radiation would get to infinity; 
in this case SEa = 0, since there is no diffusion through 
the null cone that expands outward from the initial po- 
sition R = a. It is meaningful to distinguish between 
the momentarily outgoing and ingoing radiation also in a 
curved, but asymptotically flat, spacetime. One can give 
either an operational or an analytic definition. Imagine 
a directional wave generator that sends all radiation in 
a fixed direction, when located in an almost flat region. 
(That makes sense, since it is known from analytic esti- 
mates, that the fraction of the backscattcrcd energy must 
fall off at least as C x {2m/aY, where C is of the order 
of unity - [27] , [28] and [29] . By choosing a sufficiently 
distant location one can make the diffused energy 5Ea 
arbitrarily small.) This generator, when carefully moved 
to a strongly curved region, will preserve its property 
of generating directed radiation, which can be initially 
purely outgoing (or initially purely ingoing). Alterna- 
tively, one can work out an analytic definition. Initial 
data can always be split into two parts, one "initially 
outgoing" (defined below; in the Minkowski spacetime 
that would all get to the infinity) and the other purely 
ingoing (its form is similar to the former just change 
r* —t into r* +t and some signs in the expansion - but it 
is purely ingoing in the Minkowski spacetime). We will 
show in Sec. 3, that the concept of initially outgoing 
waves is useful in bulding a nontrivial construction, and 
that fact in itself justifies this notion. 



C. Initial data for wave equations 

Let us define 

*(«,t) = *.(r--,) + 2l(^ + i2<9^. (9) 

where {i = 0, 1, 2) satisfy following relations: 
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dt^i = 3*0, 
dt^i = 3*0, 
dt^i = *o. 



dt^2 = ^i-mdt^i, 

TfX 

dt^2 = *i - -^dt^u 
*2 = 0, 



(10) 



for the polar GW, axial GW and dipolc EW, respectively. 
In the Minkowski spacetime the function * exactly solves 
Eq. (2). We assume that *j(r*, t = 0)'s vanish for R < a. 
Notice that only one of the three functions (for instance 
*o) can be freely chosen. 

We will say that initial data are purely outgoing if on 
the initial hypersurface * = * and dt'i' = 9t*. The full 
solution of Eq. (2) can be now split into the known part 
* and an unknown 6, 



* = ^1 + 5, 



(11) 



with null initial values for 6 and dtS. 6 is evolved accord- 
ing to the inhomogeneous wave equation 



{-d^ + d^,)s = vs + v, 



(12) 



where 



V ■ 




*1 *2 



+ 



for polar gravitational modes, 



i?5 



for axial GW and 



V = 6to^*i 



(13) 



(14) 



(15) 



for the electromagnetic case. The splitting (11) has 
been crucial in obtaining the analytic estimates of the 
backscattcr [27] - [30] but it appears to be advantageous 
also from the numerical point of view. 



III. EXTREMIZING THE DIFFUSION 
PARAMETER 

A. Diffusion parameter and the variational problem 

Lot us define the reprocessed radiation (RR) as 
that reaching a distant observer after the passage of the 
initial pulse; the delay is caused by multiple backscatter- 
ings. RR would be absent in the Minkowski spacetime. 
(For an example, see on Figs 8-13 the parts of waveforms 
to the right from a; = 0.) We study hereafter the RR 
generated by initially outgoing waves, in order to sepa- 
rate the genuine effects of the geometric curvature from 
those implied by artificial initial data. 



The diffusion parameter k is defined as the ratio of the 
diffused energy and the initial energy, 

.... 

Our aim in this section is to provide outgoing initial data 
that maximize k. This will be done in a class of data 
that do vanish for R < a. The intuition behind this 
is that if K is large then the fraction of the energy of 
the reprocessed radiation should be also large. That 
in turn should translate into effects like vigorous ringing 
modes or tail terms. Wc conjecture, that there exists 
a correlation between k and (defined in some way) the 
strength of QNM's. 

Expressing things in technical terms: we want to max- 
imize the nonnegative quadratic form 5Ea while keeping 
fixed the positive quadratic form E{a,A]). In numerical 
calculations this task reduces to a multi-dimensional al- 
gebraic eigenvalue problem, as we shall demonstrate. In 
the first step we choose some large Rx":^ a - the upper 
end of the initial support and maximize k in the fu- 
ture domain of dependence of (a, Ri ) with the appex at 
= R{^lLiRi)+i:i2l)^t = Obviously the 

change of i?i would change k as well, but it has bec;n es- 
tablished that above some critical value of Rx the value 
of K stabilizes. It has been found by the method of trial 
and error that the choice R\ w 150m is satisfactory. 



B. Discretization of the variational problem 

In the second step we determine a functional dis- 
crete basis {fi] {i = 1...N) on the closed interval 
[r*{a),r*{Ri)]. The dimension of that basis was usually 
250 (but tests with smaller and bigger dimensions were 
also done) a number much smaller than the number of 
points (8000) in the spatial grid; that facilitated greatly 
the numerical calculation, without loosing accuracy. The 
best results have been obtained for the basis consisting 
of the first 250 Legendre polynomials with odd indices. 

Let the expansion of the function *o (the only free 
function in the initial data set - see the remark folowing 
Eq. (10)) be 



N 



r* ^ t — r* (a) 
r*{Ri) - r*{a) 



(17) 



Then one finds from Eq. (6) that the total initial energy 
is a positive definite quadratic form, 

E{a,Ri,t = 0) = 



N 



Bij{a,Ri)CiCj, 



(18) 
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where the matrix Bij is known from the numerical cal- 
culation. 

Each element ^"0,1 = fi determines some initial val- 
ues {'^ f-)t=Q, (dt'^ f^)t=Q: they give rise to solutions 
in the domain of dependence. These solutions are lin- 
early independent (due to the uniqueness of solutions of 
Eq. (2)). Therefore the solution generated by the ini- 
tial data defined by (17) can be expressed as the linear 
combination 

N 

'^{r*,t) = Y,Ci^f.- (19) 

i=l 

Thus the energy SEa diffused through the null cone con- 
necting (a,0) with {R2,t) has the form 

N 

5E{a, R2) = R^Cfij. (20) 

Again the matrix A is obtained numerically. 

The task of maximizing the ratio of the two quadratic 
forms is equivalent to finding eigenvalues in the general- 
ized eigenvalue problem 

JV N 

Y,AjCj = Xj2BijCj, (21) 

i=i i=i 

where A is the eigenvalue and (Cj) is the corresponding 
eigenvector. There are many excellent numerical proce- 
dures for solving of the generalized problem. We choose 
one from the fast EISPACK package. This allowed us 
to find several largest eigenvalues A^'^^ the eigenvectors 
(Cj*^^) and, from Eq. (17), the corresponding functions 

¥„^\a,Ri,N) for k = 1,2.... Having ^'^^^ one finds ini- 
tial data using Eqs (9) and (10). 

As a consistency check, in number of cases the wave 
packet given by '^o{a, N) was evolved and the diffu- 
sion parameter was found directly from the definition. In 
the case of disagreement the procedure could be repeated 
with other values of numerical parameters. The disagree- 
ment was never observed for the vectors maximizing k, 
but it was found in number of cases with fourth and fifth 
eigenvectors (by convention, the eigenvectors are ordered 
according to the decreasing eigenvalue, Ai > A2 > ...). 
The parameters (A'' = 250, r*{Ri) « r*{a) + 160m, the 
size of the grid) that are reported above seem to be opti- 
mal, in the sense that the corresponding integration time 
was not too long while the accuracy was kept reasonably 
good. These values have been obtained by performing 
many series of numerical calculations. 

C. Final prepciration of extremal initial data 

These pre-prepared initial data that are maximizing 
within the chosen region (in the future dependence 



zone of data defined on (a,i?i)) undergo a process of 
extending the initial data beyond Ri . Strictly saying, we 
match a function 

/(r*) =C2 + Ciexp(-r7lO) (22) 

to each eigenvector "i!Q^\a, Ri, N). The matching is dif- 
ferentiable and the gluing point G is selected indepen- 
dently for each eigenvector. The value of G has been 
obtained as follows. Fixing a and one finds (the 
upper index is put here in order to stress the local char- 
acter of the procedure) and initial values of the locally 
extremizing solution "i!^^ . With the increase of Ri , while 
keeping a fixed, the function changes. In the limit 
one should in principle obtain the sought extremizing so- 
lution, '1/ = lim/{j^oc In the numerical practice the 
integration region must be finite. The dependence of ^0 
on Ri suggests that '^^^ w const outside some region 
of compact support. The point G is numerically deter- 
mined as being some point near the transition region. 
In our cases we obtained r*(8m) < r*{G) < r*(25m). 
Therefore the chosen ^Pq bears on an asymptotically con- 
stant value. Fig. 1 shows initial profiles of the the third 
eigenvector ^0 for a = 2.001 (GW, the polar mode). 
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FIG. 1. An eigenvector is matched to an asymptotically 
constant function (lower branch - solid line). The upper 

branch (broken line) represents the eigenvector before match- 
ing. The X-axis shows r* — r*(a) and is scaled in units of 
m. 

We would like to point that this process of matching is 
to a degree arbitrary and the obtained eigenvectors can 
be expected to be close (but not necessarily identical) to 
the extremizing eigenvectors. 

IV. NUMERICAL RESULTS 

A. Extremizing initial data and k versus a 

Figs. 2 and 3 show the distribution of the initial energy 
densities of the first and the fifth axial GW modes. As 
one can expect, the mass center is closer to the horizon 
in the case of the extremal data, while (not so obviously) 
the graph of the fifth vector suggests a larger contribution 
of high frequency radiation. 
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5 10 15 20 25 30 
FIG. 3. Axial GW. Initial energy densities for the 1st (Fig. 
2) and 5th (Fig. 3) eigenvectors. Here a = 2.01. The x-axis 
shows values of r* — r* (o) and is scaled in units of m. 

Fig. 4 demonstrates that the energy support of maxi- 
mizing initial data increases with the increase of a. The 
larger a the larger distance at which the value of the en- 
ergy stabilizes. That feature of the maximal initial data is 
counter-intuitive at the first glance, since the backscatter 
is strongest in regions with large vahics of the potential 
V (around R « 3to or r* « 0) and one would expect 
accumulation of the energy near a if a » 2m. The rea- 
son why it is not so is that the backscatter depends also 
on the frequency; a radiation accumulated at a would be 
dominated by high-frequency waves, which are weakly 
backscattered. 
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FIG. 4. Polar GW. Initial energy E{R) = E{a, 0)-E{R, 0) 
{y— axis) as a function of r* — r*(a) for a — 2.1m (solid line), 
a = 3.1m (broken line) and a = 4m (dotted line). 

The scale of the ordinate is arbitrary while the abscissa 
is in units of m. 



The main lesson that can be drawn from the foregoing 
discussion is that the extremizing initial data can occupy 
a large region that extends far away from the black hole 
horizon. 

A question arises, whether one can have modes with 
large k in the case of waves that are initially well sepa- 
rated from the horizon, i.e., when a ^ 2m. As it hap- 
pens, in order to give an answer one has to combine the 
numerical approach and an analytic insight. This is be- 
cause the numerical time is proportional (with some large 
coefficient) to (a/2m)^ and the numerics is feasible only 
when a is not too big. Fortunately, analytic estimates 
show that the diffusion parameter quickly decreases with 
the distance, at least as quickly as (2m/a)^, and becomes 
small at large a [31]. Therefore if a » 2m then no modes 
with large k can exist and it suffices to restrict the present 
analysis to a being relatively small. In this paper the nu- 
merics is done for a < 6m. 

Figs 5-7 show the dependence of k on a for five succe- 
sive eigenvectors with largest eigenvalues, in each of the 
considered wave sectors. While at a very close to 2m the 
largest eigenvalue is close to 1 in all three cases, then 
the eigenvalue fifth in the order is smaller than 0.01 for 
EW and close to 0.1 for polar GW, with the axial GW 
lying in between. The next observation that should be 
made is that with the increase of a, the largest eigen- 
value changes slower than the remaining ones and the 
the falloff of eigenvalues is quickest for EW and slowest 
for the polar GW. 
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FIG. 7. Eigenvalues k for the first five eigenvectors, dipole 
EW (Fig. 5), axial GW (Fig. 6) and polar GW (Fig. 7), in 
dependence on a. The points are connected by lines (solid, 
broken, etc.) in order to make easier the identification of 
eigenvectors. The x-axis shows the position of a and is scaled 
in units of m. 



B. On the stability of templates 

Our earlier observation that QNM can be born and can 
die [24], when observation points are moved away form 
the black hole horizon, can be rephrased as demonstra- 
tion that templates can critically depend on the distance 
of an observer from the horizon. Below we repeat that 
study and establish a lower bound on the distance of the 
observer from the horizon that is needed in order to de- 
tect a reliable wave profile. 

Fig. 8 shows that there arc many oscillations at 
R = 10m, which gradually die when the observation 
point is moved away to R= 100m (Fig. 10). One can see 
that only the first eigenvector produces some distorted 
oscillations at i? = 100m while the remaining two fail 
completely to show any ringing. 
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FIG. 9. 
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FIG. 10. Polar GW, a = 3rn. Waveforms (hi I*]) generated 
by the three strongest eigenvectors (no 1 - solid line, no 2 
- broken line and no 3 - dotted line), taken a,t R = 10m 
(Fig. 8), R = 25m (Fig. 9) and R = 100m (Fig. 10). The 
abscissa is in units of m. The scale of the ordinate is arbitrary 
and the amplitudes of each type of eigenvectors arc properly 
normalized, for the sake of clarity. The data to the right of 
a; = correspond to the purely backscattered radiation. 

Notice that while the amplitude of the surviving QNM 
seems to increase moderately, the tail (and pre-tail) part 
extends and gains in power significantly. This agrees with 
conclusions of [24]. Particularly interesting is the com- 
parison of templates shown in Fig. 9, taken at i? = 25m 
[32] and in Fig. 10, determined at i? = 100m. They are 
clearly different - the ringing phase can be much shorter 
or even disappear, while the remnants of the initial data 
(the parts of the diagrams to the left from x = 0) seen at 
R = 25m are completely different from those detected at 
R = 100m. One can conclude that the process of taking 
waveforms is unstable under the translation of the obser- 
vation point - the templates can strongly depend on the 
location of the observer. 

One can also infer from the preceding information that 
25m is too close for being the observation point and 100m 
may well be the lower bound for the observer's position. 
To this point, let us add that in many of analyzed exam- 
ples (not reported here) the waveforms did not change 
significantly above R = 100m. 
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C. strong ringing modes 

One of our aims is to find initial data that give the 
strongest possible ringing within the reprocessed radi- 
ation. The diffusion energy 5Ea bounds the energies of 
QNM, the tail (and pre-tail) term and also of the radia- 
tion falling to a black hole. While we do not have ana- 
lytic estimates of the shares of the particular contributing 
terms in 5Ea, it is obvious that configurations with large 
K have some room for robust oscillations. For that reason 
we study waves defined by the extremal initial data. 

Figs 11-13 present the radiation corresponding to 
the EW and GW initial pulses as seen by an "observer" 
situated at i? = lOOm. The x = point of the abscissa 
corresponds to the moment of time t = r* (100m) — r* (a) . 
This train of data that moves with the speed of light is 
seen earlier [t < lOOm — r*(a)) and it lies to the left from 
X = 0. To the right from a; = we have t > 100m — r* (a); 
in the absence of the backscattering there would be no 
signal at all. 
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FIG. 12. 
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FIG. 13. Templates (In 1*1) of the axial GW (Fig. 11), po- 
lar GW (Fig. 12) and EM waves (Fig. 13), the three strongest 
eigenvectors (depicted as in Figs 8 - 10). The "observer" is 
located at R = 100m and the parameter a = 2.001. The part 
to the right of a; = is the purely backscattered radiation. 

Notice that the amplitudes of the strongest ringing 
mode arc of the order of the largest amplitudes of rem- 
nants of the original radiation. This is particularly 
clearly manifested in the case of the strongest polar GW 
eigenvector. Observe also a strong deformation of the 
original signal just before t = r*(100m); initial waveform 
would be zero at x = 0, while in Figs 10-13 one can 
see a gradual build-up of a backscattered signal. Again 
the effect is strongest for the polar GW (Fig. 12, the 1st 
eigenvector) , when the backscattered part exceeds the re- 
maining signal by a factor of 10. We would also like to 
direct the attention of the reader to Fig. 10. There the 
ringing is absent for the second and third eigenvector, 
but a very strong pre-tail term is observed, comparable 
to the remainder of the main signal. 

These examples essentially confirm our conjecture that 
there exists a correlation between the diffusion factor and 
some features (strongest QNM and/or the longevity of 
the ringing phase) of the ringing. (The reservation "es- 
sentially" is caused by the fact, that the ringing belonging 
to the second axial mode on Fig. 12 is stronger than that 
of the first eigenvector; but in this case the diffusion pa- 
rameters differ only by the factor of 2.) An intuitive 
explanation with analytic flavour would be following. 
There is effectively (St^*)^ contribution to the observed 
energy flux, if an observation point is located far away 
from the horizon (the asymptotic zone, where radiation 
is dominated by the \E'o-type term and di'i' « —dr»^). 
Quasinormal modes oscillate and therefore they give a 
more significant contribution to the total backscattered 
energy than, say, tail terms. Hence small k would be pro- 
hibitive for any ringing, while strong k leaves this possi- 
bility open. This reasoning suggests also that the diffused 
energy might well be the best measure (imperfect, admit- 
tedly) of the energy of quasinormal modes generated by 
moving wave pulses. 

It was reported earlier (see, for instance. Sec. IX in 
[35]) that there exists a (sharp value of) critical width 
(suitably defined) of initial data corresponding to strong 
ringing and that both (sub- and super-) critical data gen- 
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erate much weaker ringing. While we observe a kind of a 
similar dependence, it is certainly less dramatic and no 
sharp indicator seems to be appropriate. Admittedly, we 
deal with a different situation - there is only an (initially) 
outgoing radiation, while in [35] there are both (initially) 
ingoing and outgoing components - but that probably is 
not relevant. More important can be a different shape of 
initial data - here determined by the extremization pro- 
cedure of the preceding section, while in [35] assumed to 
be gaussian. 

V. LINEAR VERSUS NONLINEAR 
DESCRIPTIONS OF THE POST-MERGER 
EVOLUTION 

A. How typical are ringing modes 

That is a basic tenet of the General Relativity, dictated 
by the beliefs in the cosmic censorship [33] and no-hair 
conjectures [34], that at some stage after plunge/merger 
a geometry generated by a pair of black holes can be 
represented as a single perturbed black hole. The per- 
turbations would be represented by gravitational waves 
and the final black hole would be either spinning (the 
Kerr black hole) or nonspinning (the Schwarzschild black 
hole), the latter in the case of the head-on collision. The 
so-called close limit approximation ( [36], [37], [38]) seems 
to assert that the linear approximation is valid since the 
formation of a common apparent horizon. Anninos et 
al. [37] give some arguments in favour of this claim that 
are supported (albeit with some reservations) by their 
analysis of the head-on collisions [39], with initial data 
of Misner type [40] . Gomez et al. [41] provide other sup- 
porting arguments in their discussion of fissioning white 
holes. An interesting new feature of a recent work by 
Husa et al. [42], which uses close approximation, is a 
weak dependence of waveforms on the collision velocity 
of two black holes. 

If this scenario is right then the naive expectation 
would be that most of the radiation is concentrated in 
the vicinity of the horizon [43] . (This is in fact observed 
- see Fig. 10 in [37], which shows that initial perturba- 
tion extends to regions very close to i? = 2m.) One can 
split these initial data into initially ingoing and outgo- 
ing parts, according to the descriptions of Sees 2.B and 
2.C. (In the example given in Fig. 10 of [37] the ingoing 
radiation remains forever inside the potential well [44].) 
The latter can be expanded in the diagonalizing basis 
defined in Sec. Ill and consisting of 250 base vectors, 
with the parameter a (that in fact specifies this basis - 
see Sec. Ill) being very close to 2m. But if a « 2m, 
then Sec. IV. A suggests that there appear a number of 
eigenvectors (from at least 2 for EW to at least 4-5 in 
the case of polar GW) with diffusion parameters being 



close to 1. The initial data for the linear phase are deter- 
mined by the preceding nonlinear evolution; if these were 
purely random, then the chance of having large n (and 
strong ringing) would be of the order of 1%. Leaving 
aside the question whether the merger phase can be re- 
garded as a random process, the least one can say is that 
the maximizing initial data of Sec. Ill should not be apri- 
ori excluded. Another argument is that the process of the 
backscatter is selective - waves longer than QNM's are 
stronger backscattered. Since the present gravitational 
wave detectors are tuned to frequences smaller (even by 
one order in the case of less massive binaries of black 
holes [36]) than QNM's characteristic for the most likely 
sources of the gravitational radiation, there are reasons 
to expect that the detected radiation will be strongly 
backscattered (even if the ringing itself would be unde- 
tectable). 

B. Taking of waveforms — some parameters 

We pointed out in Sec. IV. B that in some examples 
the templates are unstable with respect to the change of 
the "observation" point. No significant changes in the 
wave profiles have been observed beyond 100m, but the 
wave profiles at i? = 25m can be profoundly different 
from those taken at the former point. This suggests that 
the determination of templates should be done at i? ^ 
25m and that R = lOOm might the smallest required 
distance of the observer. In some cases the duration of 
the remnants of the original signal and of the comparable 
in strength backscattered part exceeds 150m (see Fig. 
12). From this one can infer that the minimal integration 
time must exceed at least 150m; it happens to be the 
longer the closer the initial pulse is located to the horizon 
of a black hole. There is no numerical calculation, in the 
case of the full nonlinear collisions, up to our knowledge, 
which satisfies both requirements. For instance in [37], 
the observations were held at ii = 25m through the time 
interval 100m. 

VI. CONCLUDING REMARKS 

The Schwarzschild spacetime can be believed to pro- 
vide a good approximation to the last phase of the col- 
lapse of two black holes - the so-called close limit ( [36] 
- [38]) bases on this idea - provided that the collision 
is (almost) head-on. At this stage the spacetime can 
be regarded as consisting of a single black hole having 
a mass m and some gravitational radiation that prop- 
agates on a Schwarzschlidean background. The initial 
data for the linear evolution of the gravitational radia- 
tion should be provided by a numerical solution of the 
preceding phases of the collapse. This task is at present 
(and presumably for some years to come) unavailable for 
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numerical relativists. In this context the existence of uni- 
versal imprints of the spacetime curvature like QNM's 
could be of relevance, but only if their amplitudes are 
strong enough. 

We invent a variational procedure that generates initial 
data corresponding to the strong backscatter. These ini- 
tial data have some features that might look as counter- 
intuitive; they have an extended support and a significant 
fraction of the wave signal energy comes from a distant 
region {R ^ 2m) . They generate strong ringing modes or 
(if the ringing is absent) robust terms preceding the tail. 
In many cases the backscattered terms and QNM's are 
much stronger than the remnants of the original signal; 
therefore the QNM's waveforms cannot be ruled out as 
objects of interest for the gravitational wave astrophysics. 

Finally, it has been shown how from the linear descrip- 
tion one can get clues as to the preparation of templates 
in the numerical analysis of the full nonliner problem of 
collapsing (head-on) black holes. 
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